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ABSTRACT 

We discuss aspects of gravitational modifications of Schrodinger dy- 
namics proposed by Diosi and Penrose. We consider first the Diosi-Penrose 
criterion for gravitationally induced state vector reduction, and compute the 
reduction time expected for a superposition of a uniform density cubical 
solid in two positions displaced by a small fraction of the cube side. We 
show that the predicted effect is much smaller than would be observable 
in the proposed Marshall et al. mirror experiment. We then consider the 
"Schrodinger-Newton" equation for an iV-particle system. We show that 
in the independent particle approximation, it differs from the usual Hartree 
approximation applied to the Newtonian potential by self-interaction terms, 
which do not have a consistent Born rule interpretation. This raises doubts 
about the use of the Schrodinger-Newton equation to calculate gravitational 
effects on molecular interference experiments. When the effects of Newto- 
nian gravitation on molecular diffraction are calculated using the standard 
many-body Schrodinger equation, no washing out of the interference pattern 
is predicted. 
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1. Introduction 

There is now considerable interest in mounting experiments to search for, and/or 
to place limits on, possible modifications of Schrodinger dynamics. We focus in this pa- 
per on conjectured gravitational modifications of the Schrodinger equation associated with 
the work of Diosi [1], Penrose [2] and their collaborators. These authors have proposed a 
gravitationally based criterion, which we refer to as the Diosi-Penrose (DP) criterion, for 
predicting when a superposition of two spatially displaced states of the same object will 
reduce to either one state or the other. In Sec. 2 we briefly review the DP criterion and its 
theoretical motivations, including the gravitationally driven stochastic equation formulated 
by Diosi [1]. In Sec. 3 we evaluate the DP effect for a uniform cube displaced by a small 
fraction of its side, and show that the predicted rate of gravitational state vector reduction 
is too small to be observed in the proposed Marshall et al. [3] mirror superposition exper- 
iment. A different, non-gravitational criterion based on displacement of the center-of-mass 
wave packet, will however be tested by the Marshall et al. proposal. 

Diosi [4] and Penrose [2] have also proposed a nonlinear equation, called the "Schrod- 
inger-Newton" (SN) equation, for including non-stochastic effects of gravitation on quantum 
evolution. In Sec. 4 we review the SN equation, give its specialization in the independent 
particle approximation, and contrast this with the standard Hartree approximation as ap- 
plied to the inter-particle Newtonian potential. We show that the two differ by a particle 
self-interaction term, which is not included in the standard Hartree approximation to Newto- 
nian dynamics, and which does not have a consistent probabilistic interpretation within the 
framework of the Born rule. Salzman and Carlip [5], motivated by searching for distinctive 
features of non-quantized gravitation, have recently argued that the SN equation implies 
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potentially observable effects in molecular diffraction experiments. In Sec. 5 we consider 
gravitational effects on molecular diffraction in standard many-body quantum theory as ap- 
plied to the inter-particle Newtonian potential, which omits the suspect self-interaction effect 
of the SN equation. We show (without invoking the Hartree approximation) that there is a 
complete decoupling of gravitational effects from the center-of-mass motion of the molecule, 
and thus no reduction in visibility of molecular interference fringes is predicted. 

2. The Diosi-Penrose (DP) criterion and Diosi's stochastic Schrodinger equation 

Diosi [1] proposed that there is a "universal gravitational white noise" , represented 
by a stochastic term <f>(r, t) in the gravitational potential (where r is the coordinate three- 
vector). Denoting the stochastic expectation by E[...j, this fluctuating part of the gravita- 
tional potential is assumed to obey 

£[0M)]=o , 

(i) 

E[<j>{r, t)0(r', *)] =TiG\r - r'^Sit - t') , 
with G the Newton gravitational constant. Including in the Schrodinger equation, Diosi 

is led to a stochastic dynamics 

ihaj,(t) = [h + J d 3 r<f>(r,t)f(r)^J if>(t) , (2a) 

with H the usual Hamiltonian and f(r) the local mass density operator. This in turn implies 
that the stochastic expectation density matrix p(t) = E[ip(t)ip(t)i] obeys the dynamical 
equation 

pit) = - ^ // yz^U^ • (26) 

Letting X denote the system coordinates, and f(r\X) the mass density at r for the 
system configuration X, Eq. (2b) implies that the off-diagonal matrix element (X\p(t)\X') 
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damps with a characteristic time Td(X,X') given by 



r d (X,X') 



G_ 

2h 



II 



d 3 rd 3 r 



Af(r\X)-f(r\X')}[f(r'\X)-f(r'\X')] 



(3) 



(Note that Eq. (12) of Diosi's paper where r d is defined contains an algebraic error, and 
should read as in Eq. (3) above, which is what one gets when one takes the off-diagonal matrix 
element of Diosi's Eq. (11). This error was noted some time ago by Anandan [6].) Although 
the density matrix evolution of Eq. (3) leads to exponential damping in time of the off- 
diagonal density matrix element (X\p(t)\X'), the stochastic Schrodinger equation of Eq. (2a) 
does not lead to state vector reduction, since an initial superposition of configurations X 
and X' does not evolve to just one of the two alternatives. However, a non-linear variant of 
Eq. (2a), constructed according to the continuous spontaneous localization scheme reviewed 
by Bassi and Ghirardi [7] and Pearle [7], does lead to state vector reduction, with the 
stochastic expectation density matrix also obeying the evolution equation of Eq. (2b). 

Penrose [2] has also proposed a role for gravitation in state vector reduction, based on 
the observation that when a macroscopic mass distribution is moved significantly, the space- 
time geometry is changed. Since standard quantum theory does not permit the description of 
coherent superpositions of states constructed on two different background geometries, Pen- 
rose argues that in a correct theory that merges spacetime geometry with quantum theory, 
such coherences must decay. He thus arrives at a criterion which states that a coherent su- 
perposition of matter density distributions p(x) and p'(x) should reduce to one or the other 
in a characteristic time r^ 1 = A/Ti, with A given by 



(In his papers, Penrose uses the notation x,y for what we have termed r,r', and his 2000 
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paper [2] giving Eq. (4) differs by a factor of Air from his 1996 paper [2]. We will follow 
the later version, and will reserve the designation x, y, z for the Cartesian components of r.) 
Apart from obvious differences in notation, and an extra numerical factor of 2, Penrose's 
criterion of Eq. (4) is the same as Diosi's criterion of Eq. (3), and we shall refer to the two 
collectively as the Diosi-Penrose (DP) criterion. 

Because Eq. (4) diverges for point particles, the effect predicted depends on the 
radius assigned to the elementary mass distributions. Moreover, the density matrix evolution 
of Eq. (2b) predicts energy non- conservation, which as discussed by Ghirardi, Grassi, and 
Rimini [8], disagrees with experimental bounds unless the point particle mass distributions 
are smeared considerably more than originally envisaged by Diosi. Rather than adding a 
smearing radius as an additional parameter of the model, we note that for any smearing 
radius greater than a typical interatomic distance of 10~ 8 cm, the mass distribution becomes 
effectively uniform. Motivated by this, we shall assume a homogeneous mass distribution in 
applying the DP criterion. 

3. Magnitude of the DP estimator in the Marshall et al. mirror experiment 

Continuing with Eq. (4), with mass distributions assumed homogeneous, let us con- 
sider the specific geometry of the Marshall et al. [3] proposal, in which a cubical mirror with 
side S = 1CT 3 cm is put into a superposition of two states displaced parallel to a side of the 
cube by d — 1CT 11 cm. Since the displacement d is a small fraction of the mirror dimension 
S, we follow Diosi [9] and Geszti [10] and expand Eq. (4) to leading, quadratic order in d. 
Writing 

p(r) = Po 9(S - x)9(x)9(S - y)9(y)9(S - z)9{z) , 

(5a) 

p'{r) = Po 9{S - x)9{x)9{S - y)9{y)9(S -z- d)9{z + d) , 
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with 9{x) the standard step function that jumps from to 1 at x = 0, we have 

p{r) - p'(r) = p o 0(S - x)9(x)9(S - y)9{y) [9{S - z)9{z) - 9{S - z - d)9(z + d)} . (5b) 
Substituting 

9(S-z- d) ~6(S -z)- d5(S - z) 

(5c) 

9{z + d) ~9{z) + dS(z) , 

we find 

p(r) - p\r) ~ d Po 6(S - x)9(x)9(S - y)9(y) [-S(z) + S(S - z)\ , (5d) 

with a similar expression with all coordinates replaced by primed coordinates. Thus Eq. (4) 
becomes 

A = Gd 2 plh , (6a) 

with 1 1 given by 

h = f f d 3 rd 3 r'9(S - x)9{x)9{S - y)9{y)9{S - x')9{x')9{S - y')9{y') 

J J (66) 

x [-5{z) + 5(S - z)\[-5{z>) + 5(S - z')][(x - x'f + (y - y'f + (z - z'f]- 1 ' 2 ■ 
Using the delta functions to eliminate the z, z' integrals, imposing the theta function 
constraints on the x,y,x',y' integrals and and scaling out the cube side S, we get finally 

A = 2Gd 2 S 3 p 2 I , (7a) 

with I the dimensionless integral given by 

I = l dX l dy lo dX 'i ( [(x ~ x'f + (y- y<n' 2 " [(x - x' f + (y - y'f + lp ) 

(76) 



The quadruple integral / can be simplified by transforming to sum and difference variables 
i] x = x — x' , a x = x + x', etc., giving the double integral form 

1 = 4 / dr * .[ - - *' iwrw " (7c) 

=2tt/3 ~ 2.0944 . 

The evaluation of the integral on the first line of Eq. (7c) was done using Mathematica®; 
as a check we also used Mathematica® to numerically evaluate the quadruple integral of 
Eq. (7b), giving the same result. 

Putting everything together, we have 

A = (A7i/3)Gd 2 S 3 p 2 , (8a) 

which with d = KT 11 cm, S = 1(T 3 cm, and p S 3 = 5 x 1(T 12 kg gives 

A =2.2 x l(T 20 uc cm" 1 , 

(86) 

r d =h/A = 1.5 x 10 9 s . 
Hence, the characteristic time for gravitational effects on the superposed cube wave function, 

according to the DP criterion, is much longer than the observation time interval of the 
Marshall et al. proposal, which is given in terms of the mirror oscillation angular frequency 
uj m by 2n/u} m = 2 x l(T 3 s. 

Thus, the Marshall et al. proposal, even it achieves the sought-for sensitivity, will 
not confront the DP proposal for state vector reduction, when interpreted using homoge- 
neous mass distributions. We emphasize at this point that the Marshall et al. paper does 
not suggest that it will test gravitationally induced reduction models (although citation of 
the Penrose papers [2] in the Marshall et al. proposal might lead readers to conclude oth- 
erwise). The mirror experiment proposal suggests a different, non-gravitational, criterion 
for state vector reduction, that superpositions reduce when an object is displaced by more 
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than the width of the center-of-mass wave packet, and this condition is met by the proposed 
experiment. The purpose of the exercise we have just gone through has been, first of all, to 
get the explicit formula for the DP criterion in the context of the mirror experiment, and 
secondly, to demonstrate that the DP criterion and the center-of-mass displacement criterion 
can make very different predictions. For completeness, we note that the mirror experiment 
may also be sensitive to other types of spontaneous localization models, if the stochasticity 
magnitude is taken large enough to give state vector reduction in latent image formation, as 
discussed in Adler [11] (which draws on earlier analyses of the mirror experiment in [12]). 

4. The "Schrodinger-Newton" (SN) equation in the independent particle 
approximation versus the Hartree approximation 

As an attempt to incorporate quantized matter into a purely classical theory of 
gravitation, M0ller [13] and Rosenfeld [14] have suggested that the source term in the classical 
Einstein equation be taken as the expectation {ip\T^ v \i/j) of the energy momentum operator 
T^ u in the quantum state \ip). As a nonrelativistic realization of this idea, Diosi [4] and 
Penrose [2] have proposed what has come to be called the "Schrodinger-Newton" equation, 
in which a quantum many-body system of N particles moves in a gravitational potential given 
by the quantum expectation of the operator Newtonian potential. Following the exposition 
of Diosi [4], the many-body equation for particles of masses mi, ...,rriN is taken as 

( N h 2 d 2 N N \ 

ihdi>(X,t)/dt= — + V rs (x r -x s )+J2^s<P(x s ,t)U(X,t) . (9a) 

\ r=l r r r,s=l s=l / 

Here V rs is a non-gravitational interaction potential, which we shall ignore for the present 
discussion, X = (xi,X2, ...,a;Ar) denotes the spatial coordinates of the N particles, and (f)(x) 
is the Newtonian gravitational potential obtained from the nonrelativistic specialization of 



the M0ller-Rosenfeld equation. In other words, is obtained by solving 

r N 

V 2 0(x,t) =4nG / d 3N X'\^(X , ,t)\ 2 J2 m uS (:i \x-x , u ) , 

J u=l 



(96) 



where X' = (x[, x' N ). Inverting Eq. (9b) and substituting into Eq. (9a), we get the 
Schrodinger-Newton equation 



N n 2 d 2 N 

ihd^(X, t)/dt = ( - + E V »(*r ~ X.) 

r=l r r r,s=l 
N 



f d ^x' mx', t )\ 2 Wo . 



(9c) 



For a single particle of mass to, this reduces to a Schrodinger equation with a nonlinear and 
nonlocal interaction term, 

ihdr/>(x, t)/dt = -^l^(x, t) - Gm 2 [ rfV ^^rW , *) • (9d) 
2to J | a; — x'\ 

Let us now specialize Eq. (9c) to the case when the non-gravitational interaction V rs 
vanishes, so that it becomes 

/ N h 2 F> 2 N r \ 

ih*Kx,t)/dt= -E^r^- G E d^x' ] p^Mx',t)\ 2 )^x,t) . (io) 

\ r=l r r u,s=l I s u\ J 

We wish to study the form taken by Eq. (10) when we make an independent particle Ansatz, 

N 

i/>(X,t) = Y[i/> r (x r ,t) , (11a) 

r=l 

with each single particle wave function ip r normalized to unity, 

J d 3 x r \i^ r (x r ,t)\ 2 = 1 . (116) 
Substituting Eq. (11a) into Eq. (10) and using Eq. (lib), and dividing by ip(X,t), we get 

N 

^%,^ s (i s ,t)=0 , (12a) 



s=l 
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with 

fi 2 d 2 N f mm 

F(x s ,t) = -ihdip s (x s ,t)/dt- — — ij s (x s ,t)-Gj2 J d3 < JpZ^, | \M<,t)\ 2 M* s ,t) ■ 

(12b) 

Since the different terms in Eq. (12a) involve independent variables x s , the usual separation 
of variables argument implies that each must be a constant, 

F(x s ,t))/ip s (x s ,t) = c s , (13a) 

with the constants c s summing to zero, 

N 

$> s = . (136) 

s=l 

However, if we introduce new single-particle wave functions ip s (x s ,t) through 

ip s (x s ,t) = exp(ic s t/h)^ s (x s ,t) , (13c) 
then by virtue of Eq. (13b), we have 

N N 

Y[ll>r(x r ,t) = Y[A(x r ,t) , (13d) 
r=l r=l 

and Eq. (13a) becomes F(x s ,t) = 0, where F(x s ,t) is obtained from F(x s ,t) of Eq. (12b) 
by replacing ip s by i\) s . We thus conclude that there is no loss of generality in taking the 
separation constants c s all as zero, and the single-particle equation as F(x s ,t) = 0, that is 

i hdMxs,t)/dt = - — —M*s,t)-Gj2 / d 3 x' mum ° \M<,t)\ 2 Mxs,t) . (14) 
2m s dxj f^J \x s -x' u \ 

Equation (14) has an almost familiar look; it has the same structure as the time- 
dependent single particle equation that one would get by treating the Newtonian inter- 
particle potential in the Hartree approximation, except that it includes a self-interaction 
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term coming from the u = s term in the summation, 




\ip s (x' s ,t)\ 2 ip s (x s ,t) 



(15a) 



Such self-interaction terms of a single particle never appear in the Hartree equation, and do 
not have a consistent interpretation within the Born rule interpretation of quantum theory. 
A term with u 7^ s in the potential energy of Eq. (14), 



has the interpretation that the gravitational potential felt by particle s at coordinate x s , as 
a result of the presence of particle u at x' u , is the Newtonian potential — Gm u m s /\x s — x' u \ 
weighted by the probability \ip u (x' u ,t)\ 2 of finding particle u at x' u . However, this interpre- 
tation does not extend to the case u = s, since when particle s is at x s , the probability 
of simultaneously finding it at x' s is zero! In terms of projectors, in the case u 7^ s we 
have that P u (x' u )P s (x s ) gives a nonzero projector for finding particle s at x s , and particle u 
simultaneously at x' u . However, in the case u = s we have P s (x' s )P s (x s ) = for all x' s 7^ x s . 

We conclude from this analysis that the Schrodinger-Newton equation does not give 
a consistent interpretation of the mutual gravitational interactions within a single system 
of particles. It can, however, be used to calculate the gravitational effect of one group 
of particles on a disjoint group of particles (say, of the sun on a planet), since then the 
problematic self- interaction terms are not present. 




(156) 
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5. Gravitational effects on molecular scattering in standard 
many-body quantum theory 



In a recent archive posting, Salzman and Carlip [5] studied the single particle case, 
Eq. (9d), of the SN equation, and based on this suggested that there may be significant 
nonlinear gravitational effects in potentially observable situations, such as molecular inter- 
ferometry experiments. However, the single particle case of the SN equation consists of a 
self-interaction term which, as we observed in the preceding section, does not appear in the 
standard Hartree approximation, and which does not have a Born rule interpretation. This 
makes it problematic, we believe, to apply the SN equation to the mutual gravitational in- 
teractions within a system of atoms, as needed, for example, to discuss gravitational effects 
in molecular diffraction. 

There is a standard way of treating gravitational effects on large molecules, within 
conventional many-body theory (without use of the Hartree approximation), which leads to 
a different conclusion from that reached in [5]. One simply includes in the interaction term 

N 

^2V rs (x r -x s ) (16a) 

r,s=l 

of Eq. (9a) a Newtonian gravitational potential term 

r^ts 

in analogy with the usual treatment of the inter-particle Coulomb potential. Since both 
Eq. (16a) (which includes the Coulomb force terms) and Eq. (16b) (which is the gravitational 
perturbation) depend only on the relative coordinates x r — x s , they do not appear in the 
equation for the center-of-mass wave function of the molecule. The center of mass will thus 
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obey a free-particle Schrodinger equation, subject to external potentials (such as diffraction 
gratings) acting on the molecule. Therefore, one expects no significant effect on the molecular 
interference pattern from the mutual gravitational interactions of the molecular constituents. 
Such gravitational perturbations will very slightly change the shape and energy levels of the 
molecule, but will not exert an influence on its center-of-mass motion, other than (when 
relativistic effects are taken into account) through their small modification of the mass of 
the molecule. 

We conclude with a question that suggests further work. As noted above, the SN 
equation arises from applying the M0ller-Rosenfeld semiclassical equation to the Newtonian 
interaction of a many-particle system. Do the problems that we have encountered indicate 
that a semiclassical approach to gravitation is inconsistent, and hence that gravity must be 
quantized [15]? Or do they only indicate that a modification of the M0ller-Rosenfeld and 
SN approach should be sought, which will make possible a consistent semiclassical theory of 
gravitational effects? 
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